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ABSTRACT.  A  method  for  recomatnciMQ  the  complex  index  of  refraction  of  a  bounded  inho¬ 
mogeneous  object  from  measured  scattered  field  data  ts  presented.  Some  numerical  examples  are 
given  indicating  the  limits  on  the  contrasts  vhich  can  be  reconstructed. 


1.  INTRODUCTION 

Assume  that  an  inhomogeneous  obstacle  D  is  irradiated  successively  by  a  number  of  known  inci¬ 
dent  fields  i  =  1,  •  •  -  ,/.  For  each  excitation,  the  direct  scattering  problem  may  be  reformu¬ 
lated  as  the  domain  integral  equation 

=  t«i(p)  -  GdX«i(p).  p  e  D  ,  (1) 

where 

GDX*^i{p)  =  J^G(pyq)x{q)Mq)dvq,  p^D.  (2) 

Here,  tii  is  the  total  field,  k  is  the  wavenumber,  x  “  complex  contrast  (x  =  «*-!»  where  n 
is  the  index  of  refraction),  G(p,q)  is  the  free-tpace  Green’s  function  and  p  and  q  are  position 
vectors.  Go  is  an  operator  mapping  L^(D)  (square  integrable  functions  in  D)  into  itself.  If  5 
is  a  surface  enclosing  D  then  the  scattered  field  on  S,  is  pven  by  GsX^i  ^here  Gs  is  the 
«*in*  operator  defined  in  Eq.  (2),  except  the  field  point  p  now  lies  on  5.  Hence  Gs  i*  an  operator 
mupping  L^{D)  into  L^{S).  We  assume  that  is  measured  on  5  and  denote  by  fi{p),peS,  the 
measured  data  for  each  excitation  t,  s  =  1,  •  •  > ,  /.  The  profile  inversion  problem  is  that  of  finding 
X  for  given  A.  or  solving  the  equation 

Gsxm(p)  =  P^s,  (3) 

for  X  subject  to  the  additional  condition  that  n,  and  x  satisfy  Eq.  (1)  in  D.  We  will  seek  Ui  and 
X  simultaneously  to  minimise  the  error  on  D  in  satisfying  Eq.  (1)  and  the  error  on  5  in 
satisfying  Eq.  (3). 


2.  THE  INVERSION  ALGORITHM 

Here  we  propose  an  iterative  inversion  algorithm  which  incorporates  the  ideas  of  successive  over¬ 
relaxation  as  well  as  the  coqjngate  gradient  method.  Specifically  we  propose  the  iterative  con¬ 
struction  of  sequences  and  {xn}  ••  foUows: 

=  «;-“"  xo  =  3r“-', 

*i»n  “  n»  Xn  =  Xn-I  +  0n^  i 

ri,n  ^  ”  GsXn^,n  9  {*) 

where  On  0n  general  complex  constants  which  are  chosen  at  each  step  to  minimise 

1=1  isl  \3=1  /  M=1  / 


and  the  subscripts  D  and  5  on  the  norm  ||  •  ||  and  inner  product  (•,•)  in  I?  indicate  the  domain 
of  integration.  The  minimization  of  the  quantity  Fn  of  Eq.  (5)  leads  tc  a  nonlinear  problem  for 
the  coefficients  and  0n  at  each  step,  which  we  solve  using  a  conjugate  gradient  method.  The 
starting  value  for  q„  is  obtained  by  taking  =  0  and  minimizing  Fn,  while  the  starting  value  for 
/3„  is  found  by  setting  Qn  =  0  and  again  minimizing  F„. 


5.  INITIAL  GUESS  AND  CORRECTION  DIRECTIONS 

In  our  previous  treatment  of  this  problem  [1]  we  chose  xo  =  0,  while  the  update  direction  for 
the  field  was  directly  adapted  from  the  successive  over- relaxation  method  for  solving  the  direct 
problem  with  known  contrast  to  be  t;,,„  =  and  the  update  direction  for  the  contrast  was 

chosen  to  be  the  gradient  of  the  error  in  the  measured  data  at  the  previous,  (n-l)st,  step.  In  the 
present  work  we  refine  these  choices  considerably. 

We  obtain  the  initial  guess  by  finding  the  constant  contrast  and  associated  fields  uj”***“^ 

which  minimize  the  fimctional  Fn.  Specifically,  we  proceed  as  foUows.  Define  the  normalized 
change  in  the  field  by 

=  (i:  II?,)  *  (e  iid)  *  •  (6) 


We  set  an  arbitrary  switching  criterion,  e,  and  run  the  algorithm  of  Eq.  (4)  with 

=  uj”^,  dn  =  1  and  until  <  c,  then  switch  the  definition  of  to 

V,,n=9ln  +  ‘l>,.n-l,  7;^  =  (^5.%  /  (E  IK»-1 1^^^^  . 


mfia/  _ 


with  the  gradient 


where  the  overbar  denotes  complex  conjugate  and  Gs  i»  a  map  from  L*(5)  to  L^{D).  The  choice 
of  the  direction  v,,n  In  Eqs.  (7)  -  (8)  is  the  Polak-Ribim  conjugate  gradient  direction  assiuning  the 
contrast  does  not  change.  Continue  this  algorithm  until  we  again  achieve  £„  <  e.  The  resulting 
values  are  taken  as  and  x*”***^^ 

With  these  initial  choices  we  run  the  algorithm  of  Eq.  (4)  with  Vi^n  ^  In  Eqs.  (7)  •  (8)  and  dn 
is  taken  in  one  of  the  two  ways:  if  £n>i  ^  ^  then  is  taken  to  be  the  gradient  direction 

gi  -  -WD^Ui^n-^iGori^n-X  +  ti;s5Z^.n-|GsPi,n-l  ,  (9) 

i=rl  1=1 

whereas  if  £„^|  <  e,  we  use  the  Polak-RJbiere  conjugate  gradient  direction 

d.  =  gt  +  lidn., ,  yi  =  ({gi,gi-gtr)D)  /  (||tf^-,||?,)  •  (10) 

Continue  the  iteration  until  either  F„  meets  a  preset  error  criterion  or  ceases  to  change. 


4.  NUMERICAL  EXAMPLES 

The  inversion  method  is  illustrated  in  a  particular  case  in  2-D  scattering  by  a  square  cylinder  of 
dimension  dxd  with  sinusoidal  varying  profile,  x  =  sin(ir£/A)sin(xy/A)  for  0  <  c, y  <  3A,  so  that 
kd\xtnai\  -  Results  are  shown  for  twenty  and  thirty  equally  spaced  measurement  stations 
distributed  on  a  circle  of  radius  3A  containing  the  cylinder  with  each  station  serving  successively 
as  the  source  and  all  stations  serving  as  receivers,  I  =  20,30.  The  cylinder  is  discretized  into 
29  X  29  subsquares.  The  original  profile  is  illustrated  in  Fig.  la.  Using  c  =  0.01  the  reconstructed 
profile  is  shown  in  Figs,  lb  and  Ic,  employing  twenty  and  thirty  stations,  respectively. 


Figure  1.  The  original  profile  (a)  and  the  reconstructed  profiles  for  20  measurement  stations  (b) 

and  for  30  measurement  stations  (c). 


The  reconstructions  shown  in  Figs,  lb  and  Ic  are  the  results  after  128  iterations  of  which  56  were 
required  to  obtain  the  initial  guess  for  20  stations.  For  the  case  of  30  stations  only  38  of  the  128 
iterations  were  needed  to  obtain  the  initial  guess.  The  values  of  the  functional,  Eq.  (5),  which 
was  to  be  minimized  were  F\2h  =  0.013  for  /  =  20  and  Fi28  =  0.002  for  I  =  30. 

5.  CONCLUSIONS 

An  iterative  method  for  complex  profile  construction  has  been  described  and  tested.  The  method 
combines  the  features  of  successive  over-relaxation,  gradient  and  conjugate  gradient  methods  to 
minimize  a  functional  consisting  of  normalized  errors  in  satisfying  the  field  equation  and  the 
error  in  matching  the  measured  data.  The  field  equation  serves  as  the  regularizer  for  the  ill- 
posed  problem  finding  a  function  in  D  to  minimize  the  error  in  solving  Eq.  (3).  The  nonlinear 
optimization  problem  is  not  linearized,  however,  the  two  components  of  the  functional  in  Eq.  (5) 
are  treated  somewhat  separately.  The  algorithm  was  constructed  to  delay  large  changes  in  the 
contrast  until  the  field  was  somewhat  stable.  This  was  the  motivation  for  the  separate  treatment  of 
the  initial  guess  as  well  sts  the  subsequent  switching  in  the  algorithm  based  on  the  size  of  €n»  The 
numerical  restilts  presented  here  as  well  as  additional  experiments  indicate  that  the  algorithm 
successfully  reconstructs  complex  contrasts  for  4d|Xma*l  £  6^*  To  achieve  reconstructions  for 
large  values  of  Xmax  ww  frequency  measurements  will  not  suffice  to  give  reasonable  resolution. 
Future  work  is  directed  toward  extending  the  method  to  include  measurements  at  more  than  one 
frequency  to  accomodate  larger  contrasts. 
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